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Abstract. In this article we review the thermodynamics of liquids in the framework 
of the inherent structure formalism. We then present calculations of the distribution 
of the basins in the potential energy of a binary Lennard-Jones mixture as a function 
of temperature. The comparison between the numerical data and the theoretical 
formalism allows us to evaluate the degeneracy of the inherent structures in a bulk 
system and to estimate the energy of the lowest energy disordered state (the Kauzmann 
energy). We find that, around the mode-coupling temperature, the partition function 
of the liquid is approximated well by the product of two loosely coupled partition 
functions, one depending on the inherent structures quantities (depth of the basins 
and their degeneracy) and one describing the free energy of the liquid constrained in 
one typical basin. 
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1. Introduction 

The potential energy of a system composed by interacting atoms is a complicated 
surface in a 3iV-dimension space. The motion of the system can be thought of as a 
trajectory over such a potential energy surface (PES). At different temperatures, the 
system explores different parts of the PES, according to the Boltzmann weight. The 
idea of focusing on the PES for understanding the physics of glass forming liquids can be 
traced back to the seminal (but talkative) paper of Goldstein He suggested that the 
dynamics of deeply supercooled liquids can be described in terms of a diffusive process 
of the system between different PES basins. At low temperatures the dynamics slows 
down since the liquid becomes trapped for a long time in a basin. 

The concept of a basin in configuration space was formalized by Stillinger and 
Weber 0, who introduced a recipe, very well suited for numerical analysis, for 
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partitioning the PES into disjoint basins. The set of points in configuration space 
connected to the same local minimum via a steepest-descent trajectory defines uniquely 
the basin associated with this local minimum. Stillinger and Weber named the structure 
of the system in the minimum inherent-structure (IS) and the value of the PES at the 
minimum inherent-structure energy (e/5). 

The increased computational facilities have significantly improved the early efforts 
of studying the PES. Nowadays, an exhaustive search for all IS has been performed for 
clusters and complete maps of the inherent structure energies are available for several 
potential modelsf^. For clusters, as well as for small proteinsQ], the connectivity 
between all IS has also been evaluated, to provide a very informative map both of 
the thermodynamics as well as of the dynamics in these small systems^]. Small size 
systems, composed by 30 to 50 atoms with periodic boundary conditions, have also been 
studied in detail recently and almost exhaustive enumerations of all IS energies are now 
available]^, ||, |^. 

In this article we review a theoretical framework in which the IS results can be 
interpreted in a convenient way (Sec. |]) and discuss the approximations requested 
for a factorization of the partition function in two functions, one describing the 
thermodynamics of the IS sub-system and one describing the thermodynamic of the 
exploration of one representative basin (Sec.|^) In the following two sections we present 
calculations of the temperature (Sec. |^) and IS energy (Sec. |) dependence of the 
configurational entropy for a bulk system. Such calculations allow to quantify the 
properties of the PES for model systems and to probe the validity of the factorization 
approximation. New information on the equilibrium and aging dynamics of supercooled 
simple liquids is provided by the presented results. An account of the results has been 
reported in Ref.[|Tn 



2. Theory 

This section reviews the thermodynamic formulation proposed by Stillinger and 
Weber 0, focusing on the concept of basins in configuration space. 

The partition function of a system composed by N identical atoms of mass m, after 
the integration over the momentum variables is 

Zn = A-=^^ I exp(-/3y(r^))rfr^ (1) 

where A = (3 /2T[m. The integral over the configuration space can be separated 
into a sum over all distinct basins 

Zn = \-'''Y.^M-P^o.) I exp(-/3A,(r^))dr^ (2) 

where Ra is the set of points composing the basin a, $q, is the potential energy of the 
minimum a and the non-negative quantity AQ,(r'^) measures the potential energy at a 
point belonging to the basin a relative to the minimum. By classifying the minima 
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according to their IS energy e/^, the sum over the basins can be separated in a sum over 
all possible values of e/5 and a sum over all basins a' with the same e/5 value. 



.^ = A Y.^M-M E / exp(-/3A„,(r")) t/r^ (3) 



Following Stillinger and Weber, we introduce an IS density of states ^(e^.,), which counts 
the number of distinct basins with IS energy between e^^ and e^^ + de^g and define a 
basin free energy e^^) as the average e/5 basin free energy, according to 

- '° ( ^1^, ) 

If all basins with the same e/5 energy have the same statistical properties then e^^) 
coincides with the free energy of a system constrained to sample only one basin and 
which does not know of the existence of the other equivalent de^gVL^ej^) basins. Zjsi can 
be expressed in terms of PES quantities, as 

Zn = J dejsVt{ejs) exp{-/3eis - /3f{/3, e/5)) (5) 

Performing a maximum integral evaluation of the partition function, the free energy F 
of the system can be expressed in the thermodynamic limit as 

F = e}s-TS,onfie*js) + fiP,e}s) (6) 

where e*(T) is the e/5 value which maximizes the integrand and Sconfi^is) = 
/cb ln((5e/5f2(e/5)). If we separate now f{P,e) in its energetic Ub and entropic Sb 
contribution, we immediately notice that the entropy associated with the basin 
degeneracy, SconfiT), can be calculated as the difference between the system entropy 
and Sb, the entropy of the system constrained in an IS e*. In the present formalism, if 
one excludes from the sum in the partition function the crystalline IS, then one can 
identify F as the fluid free energy for all T. 

The choice of separating the liquid free energy in a sum of two inter-related 
contributions (via the e/5 dependence of /) has been often used in the past, for example 
in estimating the configurational entropy from available experimental data|jll|. In this 
case, the basin entropy is identified with the entropy of the stable crystal at the same 
thermodynamics point. Such an identification is based on the idea that the vibrational 
properties of a system constrained in a deep basin are similar to the properties of the 
close crystalline structure. More recently, the consequences of such separation for several 



thermodynamic quantities have been explicitly worked out 1 12]. For recent related work 
see also O, HI 



3. Low T approximation 

There are two interesting cases which may help understanding the low T dynamics of 
liquids. These cases are connected to specific forms of e/5). The first describes the 
case where e/5) ~ f{P), i-e. there is no T-dependence through e/5. In this case 
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basins are characterized by approximately the same shape in configuration space, an 
hypothesis which can be tested by studying the cjs dependence of the density of states. 
In this approximation Eq. (|^) factorizes in 

Zn ~ • (7) 



where 



and 



Z^^ = I n{e,,)exp{-/3eis)deis (8) 



Z' = exp(-/3/(/3)) (9) 

In the range of T where this approximation holds, the system can be tough of as two 
weakly coupled subsystem: the JS'-subsystem, which has now been transformed in a 
system with levels labeled by the cjs value with degeneracy fl{eis), and the basin 
subsystem which describe the motion in the characteristic basin. The coupling between 
the two subsystem, which of course allows for the equilibration process between the two 
subsystems, is due to the weak T dependence of cjs, which is neglected in the present 
approximation. 

The second case is the case when (3f{[3, ejs) ~ g{l3) + h{eis), i-e. when the T and 
Cjs dependence are not mixed. This case is realized for example in the case where at 
different T the system populate basins which are always harmonic, but different in their 
density of states[|lj|. In this second factorization of Z^ as in Eq. is also 

possible by redefining the density of states to include the basin volume in configuration 
space, as fl{ejg)e^^'''^'^^ . 



Analysis of computer simulation data allows to look for the existence of a T 
range where the factorization approximation holds. Indeed, the probability density 
of extracting from a system in thermal equilibrium at temperature T a configuration 
belonging to a basin with IS energy e/5 is 

Pics, T) = ^ie,,)eM-f3eis-mf3,e,s)) ^^^^ 

Zn{T) 

If the factorization approximation holds, then the only cjs dependence in the rhs of 

ln[P(e,„T)<5e,,] +/3e,, = S,onf{eis)/kB - HZnIP)] - (3f{(3,ejs) (11) 

is contained in Sconf- This imply that curves for different T can be superimposed after 
a shift of a T-dependent quantity. The resulting cjs master curve is, except for an 
unknown constant, the ejs configurational entropy. 



4. System 

We have studied the well-known 80-20 Lennard- Jones A — B binary mixture (BMLJ), 
composed by 1000 atoms in a volume Vo = (9.4)^, corresponding to a reduced density 
of 1.2039. Units of length, energy and are defined by the a and e of the A — A Lennard 
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Jones interaction potential, and the unit of mass by the mass of atom A. The pair 
potential is defined in Ref.|T5|. The equilibrium and out-of-equilibrium slow dynamics 
has been studied extensively. The critical temperature of Mode Coupling Theory for 
this system is 0.435|T^ . 



New simulations, covering the range 0.446 < T < 5.0, have been performed 
in the canonical ensemble by coupling the system to a Nose '-Hoover thermostat . 



From simulations lasting more than 60M steps, we have extracted 1000 equally spaced 
configurations and we have calculated for each of them the corresponding IS. 

5. Temperature Dependence of the configurational entropy 

We have estimated the T-dependence configurational entropy for the BML J as difference 
of the liquid entropy and of the basin entropy, as discussed in section 0. An independent 
(and consistent) estimate of the same quantity in the same system has been performed 
with a similar procedure by Coluzzi et al. [T^. For related work see also Refs. |TP[ 
The entropy of the liquid has been calculated via thermodynamic integration 
starting from the ideal gas binary mixture reference point, along the T = 5.0 isotherm, 
up to the studied p = 1.2 density. In the following we call C the state point 
T = 5.0, p = 1.2. The entropy of the liquid in C can be written as 

S{C) = + ^ + £ ^ (12) 

where 

S^^-^-'^-(T^P) = In(^) - ^ In(^) + ^ In(^) - ln(^)(13) 

Pex is the excess pressure over the ideal gas value, U is the potential energy, and e is the 
energy per particle. Fig. ^ shows the excess pressure as a function of the volume 
calculated from twenty-six independent molecular djTiamics simulations. At large 
volumes, the calculated excess pressure coincides with the first correction to the ideal 
gas law which can be analytically calculated from the first virial coefficient of the binary 
mixture, B2{T), which in the case of our system is equal to B2{T = 5.0) = 0.53622. To 
decrease the numerical integration error we analytically calculate the contribution to 
the integral arising from the first virial corrections and integrate over the volume only 
Pecc - B2{T)kBT{N/Vf. As a result, we obtain S{C)/kB = 8061.70 

The entropy at any T along the studied isochoric path can then be calculated as 

^(T, p = 1.2) = SiC) + f ^^dr (14) 

Jt=5.0 1 

where Cy(T) = dU{T)/dT + 3/2NkB is calculated from the T-dependence of the system 
average potential energy U obtained from the simulations. We find that, in agreement 
with recent theoretical predictions [^], the T dependence of U along the studied isochore 
is very well described by the law U (T) ~ T^/^ (See Fig.^, which produces a contribution 
to the liquid entropy varying as T~^/^. The use of the T~^/^ law provides a reliable 
extrapolation of Suquid below the lowest studied T. 
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Figure 1. Excess pressure as a function of the volume at T = 5.0. The dashed hne is 
the analytically calculated first virial correction to the pressure. 
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Figure 2. Potential energy as a function of T^/s ^-^c BMLJ system. The inset 
show the same data in liner scale. The continuous line is the &t to U = a + bT^^^ 

To estimate the basin entropy, we assume that at the lowest studied T, the 
unknown e) can be approximated by the harmonic free energy of a disordered 
system characterized by the eigenfrequencies spectrum calculated from the IS at the 
corresponding T. In this approximation, the difference between the entropy of the liquid 
and the entropy of the harmonic disordered solid coincides with the configurational 
entropy. We evaluate the entropy of the disordered solid in the harmonic approximation 
as 

3N-3 

Sdisordered-soHd{T, V) = ^ ^ — ln(/3^a;j) (15) 

where cuj is the frequency of the j-th normal mode. The e/^-dependence in 
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S disordered- soUdiT , V) eiiteis via the 6/5- dependence of the density of states. Consistently 
with the estimate of the cjs dependence of f{ejs,T) discussed below, we find (see Fig. 
^left) that the T dependence of the density of states accounts for only a few percent of 
the T-change in S disordered- solid at low T. 

The T dependence of the evaluated liquid and disordered solid entropies is reported 
in Fig.^left. The T dependence of the configurational entropy (difference between 
Siiquid and S disordered- solid) IS reported in H-middle. We note that, if the extrapolations 
are reliable, the configurational entropy vanishes at T = 0.297 ± 0.01, which defines 
the Kauzmann temperature T/^[^ for the studied binary mixture, in agreement with 
the findings of Coluzzi et alfl^, |T^. Note that the configuration entropy around 
Tmct = 0.435 is halfway between Tk and the high T value, suggesting that the ordering 
process in configuration space at the lowest temperature which we have been able to 
equilibrate is far from being complete. Of course, the present data do not furnish a 
full proof of the existence of a finite T at which Sconf goes to zero, being based on a 
large (but apparently reliable, see Fig.|^) extrapolation in T. We also note that the 
ratio between Tk and Tmct support the view that the studied system has intermediate 
fragility character, as recently predicted by Angell and coworkers on the basis of a 



comparison between experimental results and numerical data for the same system[23 



6. IS-energy dependence of the configurational entropy 

In this section we show that in the BMLJ case, for T < 0.8, the factorization 
approximation discussed in Sec. ^ is indeed satisfied. The possibility of separating 
the IS subsystem thermodynamics from the basin thermodynamics allow us to calculate 
the e/5 dependence of the configurational entropy and thus to estimate the number of 
basins in configuration space with the same ejs energy. 

To test the validity of the factorization approximation, we evaluate the Ihs of Eq. 
PI , i.e. we calculate the eis dependence of ln(P(e/s', T)) + eis/T. As discussed in Sec. |, 
if f{P, cjs) has only a weak dependence on cjs, then it must be possible to superimpose 
curves at different temperatures which overlap in cjs- Then, the resulting e/s-dependent 
curve is, except for an unknown constant, Sconfi^js), in the cjs range accessed at the 
studied T. 

This procedure is displayed in Fig.^. We note that while below T = 0.8 curves for 
different T lie on the same master curve, above T = 0.8, curves for different T have 
different cjs dependence, highlighting the progressive e/^-dependence of f{P,eis)- 

The data presented in Fig. ^ are particularly relevant. They show that, below 
T = 0.8, the IS can be treated as a system of levels characterized by an energy value 
eis and an associated degeneracy fl{eis)- Thus, for the ejs subsystem it is possible to 
use the standard thermodynamics relations to evaluate the T-dependence of the average 
energy and entropy. In this respect, the T dependent configurational entropy (but only 
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Figure 3. Left: Liquid (circles) and disordered-solid (squares) entropies as a function 
of T. The diamonds show the T dependence of the disordered sohd entropy, once 
the exphcit T dependence is subtracted, to highhght the weak residual T dependence 
due to the T-dependence of the density of states. Such weak residual T-dependence 
has been extrapolated to lower T via a quadratic fit and used to provide an analytic 
expression for the T dependence of the disordered solid entropy. Middle: T-dependence 
of the configurational entropy. Right: T-dependence of the IS energy for the BMLJ 
system as determined from the simulation (circles) and from Eq. (|l6|) (solid line). 



below T = 0.8) can be evaluated as 

^^^^ = ^ e,s{T) = ejsiTK) + f TdS^onfiT) (16) 

By integrating the configurational entropy from Tx upward, it is possible to calculate the 
the T dependence of ejs- The unknown integration constant eisiTx) can be calculated 
by comparing the obtained expression with the e75(T)-dependence calculated directly 
from the simulation in the region T < 0.8 (see also Ref. |p4|| ) The present analysis (see 
Fig. |-right) predicts eisiTx) = -7.82 ± 0.01. 

From S{T), evaluated in the previous Section, and from cjsiT), evaluated 
according to Eq. (|T^), it is possible to eliminate the T dependence and to calculate 
the ejs dependence of the configurational entropy in an absolute scale, which can be 
compared with the one calculated independently — but with an unknown constant — 
via the superposition of the different \ia{P{eis,T)) + ejs/T curves. Such a comparison 
is done in Fig. ^. The agreement between the two set of measurements confirms the 
validity of the analysis presented in this article and the quality of the factorization 
approximation. 
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Figure 4. Top: Distributions P{eis,T) of the IS energy (per atom) 
for different equilibrium temperatures T. From left to right: T = 

0.446, 0.466, 0.5, 0.55, 0.6, 0.8, 1.0, 2.0, 4.0. Middle: ln[P(e7s, T)]+^eis, for six different 
equilibrium temperatures T. Same symbols as in the top panel. Bottom: Data in the 
middle panel (plus data corresponding to other T) are displayed shifted to maximize the 
overlap between curves with different T and the overlap with Sconf{s-is) (in absolute 
units), calculated as discussed in the text and shown here as full line. The curves 
which do not lay on the continuous line correspond to T = 5.0,4.0,2.0,1.5,1.0,0.8, 
from bottom to top. 



Before concluding this section, we note that an estimate of the eis dependence 
of the configurational entropy, based on the analysis of experimental data have been 
presented in Ref . [p5|] . Analysis of the configurational entropy as a function of internal 
system parameter's (which conceptually are equivalent to the eis choice adopted in the 
present work) have been reported in Ref. |2^, ^ 
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7. Conclusions 

The data and the analysis reported in this article offer a detailed thermodynamic 
description of the supercooling state, based on the formalism proposed by Stillinger 
and Weber. In particular we have presented a quantitative evaluation of the degeneracy 
of the inherent structures (which before was only calculated for systems composed by 
less than 50 atoms 0, H) for a bulk system. We consider particularly relevant the 
presented evidence that in the supercooled states (below T = 0.8 for the studied 
system) the thermodynamics of the inherent structures almost completely decouples 
from the "vibrational" thermodynamics (i.e. from the process of exploration of the IS 
basin). It is particularly important to notice that a thermodynamics approach for the 
inherent structures subsystem becomes possible in supercooled states. The description of 
supercooled liquids as composed by two weakly coupled subsystems — the IS subsystem 
and the "vibrational" basin subsystem — offers stimulating ideas both for a microscopic 
understanding of the out-of-equilibrium thermodynamics recently proposed (since if 



the factorization were exact, one could think of keeping the two subsystem coupled to 
two different temperatures) and the aging processes ||2^, as well as for the still missing 
theoretical quantitative description of the slow dynamics below the MCT temperature. 
A first step in the direction of estimating the temperature at which the configurational 
subsystem is in quasi-equilibrium during an aging process has been recently reported 



Finally, we stress that the description we have presented refers to a constant volume 
system. In this respect, it is based on one internal parameter only (in the language of 
Davies and Jones[^), which we have identified with cjs- In a full treatment, at least one 
other internal parameter would be necessary, to discriminate between basins with the 
same cjs but different volume. We plan to further test the validity of such one-internal 
parameter description for isochoric cooling. 
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